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We study how Co- and Ni-substitution affect the electronic structure of the iron-based super- 
conductor, LaFeAsO. We perform ab initio supercell calculations and unfold the first Brillouin 
zone (BZ) to calculate the spectral function in the BZ for the normal cell. The charge density 
distribution in real space shows that doped extra electrons are trapped around Co (Ni) atom. 
This seems to suggest that carriers are not doped by Co(Ni)-substitution. However, the present 
momentum-space analysis indicates that the Fermi-surface volume indeed expands by substi- 
tutions, which can be well described by the rigid-band shift approximation. By taking into 
account this effective doping, we discuss whether the sign-reversing s-wave (s±-wave) scenario 
is compatible with experiments. 
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Since the discovery of superconductivity in F-doped 
LaFeAsO, iron-based superconductors are now attract- 
ing broad interests due to their high-transition tempera- 
ture (Tc) up to 56 K.^-* It has been experimentally found 
that we can enhance or even induce the high Tc supercon- 
ductivity by several ways such as chemical substitution, 
(intentional) introduction of defects, application of exter- 
nal pressure, and so on. Among them, transition-metal 
substitution is of great interest and has been intensively 
studied both theoretically and experimentally, since it 
has two important effects from which we may get crucial 
hints to clarify the pairing mechanism: One is carrier 
doping effect and the other is pair breaking effect. 

Theoretically, it has been shown that the local po- 
tential of transition metal ions has a systematic chem- 
ical trend, and the position of the Fermi level (Ep) in 
the density of state (DOS) exhibits a monotonic shift 
as a function of the concentration of extra electrons do- 
nated by the transition- metal substitutes.^' On the 
other hand, the charge density distribution p(r) of ex- 
tra electrons are localized around the impurity site,^^ 
which seemingly suggests that transition-metal substi- 
tution does not dope mobile carriers into the system. 
While these two observations seem to contradict with 
each other, it should be noted that the structure of 
p(r) does not necessarily give conclusive information on 
whether the extra electrons are localized or not. Due to 
the local charge neutrality condition, p{r) always has a 
peak structure around the impurity site, whose ampli- 
tude is determined by the Friedel sum rule. Experimen- 
tally, systematic analysis of transport properties suggests 
that carriers are indeed doped into the system.^' Co-, 
Rh-, Ni- and Pd-doping dependence of Tc were also in- 
vestigated systematically, and it was found that Tc seems 
to depend only on the amount of extra electrons. This 
behavior can be understood if the charge-doping effect 
is appreciable, but the pair-breaking effect is negligibly 



small. 

A standard way to study impurity effects by first- 
principles calculation is to introduce a large supercell 
and replace one of the atoms in the cell with a dopant 
atom. Although this approach can mimic a randomly 
doped system when the size of supercell is sufficiently 
large,*) it also has a drawback; the first Brillouin zone 
(BZ) for the normal cell are folded into the small BZ, 
where many bands are heavily entangled. Regarding this 
problem, recently, a technique of BZ-unfolding was devel- 
oped, which enable us to transform the band structure 
in the original small BZ for the supercell to that in the 
BZ for the normal cell. 

In the present study, using this method, we examine 
the effect of Co- and Ni-doping on the electronic struc- 
ture. By visualizing the evolution of the Fermi-surface 
volume by transition-metal substitution, we show that 
the system is effectively doped, although the charge den- 
sity of extra electrons are localized around the transition 
metal impurity ions. We estimate the amount of doped 
carriers and show that this doping is well described by 
a rigid-band shift of the Fermi level. Finally, we discuss 
the impurity effect on the superconductivity. 

Let us move on to the detail of the method em- 
ployed in the present study. First, we performed ab 
initio density-functional calculations with generalized- 
gradient-approximation (GGA),^°) using Tokyo Ab ini- 
tio Program Package^^^ with a plane wave basis set and 
norm-conserving pseudopotentials.^^^ The cutoff energy 
for the wave functions and charge density were set to 
be 64 Ry and 256 Ry, respectively. We replaced one Fe 
atom in 2x1x1, 2x2x1 and 3x3x1 supercell with Co or 
Ni.*) Since the normal cell contains two Fe atoms, these 
setups correspond to 25 %, 12.5 % and 5.6 % doping, 
respectively. While a fc-mesh of 5 x 5 x 5 was used for the 
former two supercells, a 3x3x3 fc-mesh was used for the 
latter supercell. 
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Next, we constructed maximally localized Wannier 
functions (MLWrs)^^-* for the Fe 3d bands around the en- 
ergy range of —3 eV to 3 eV. We then represent the Kohn- 
Sham Hamiltonian "H^^ using the resulting supercell 
Wannier orbitals as H^%,{Ii,R') = (7VR|H^^|iV'R'), 
where |iVR) denotes the A^-th Wannier orbital sitting at 
the lattice R. Note that ■H^%,{'R, R') = H^%,(R - R') 
due to the translational symmetry. The advantage of the 
MLWF method is that we can directly estimate micro- 
scopic quantities in real space, associated with the im- 
purity effect. '^^ 

Table I summarizes estimated onsite impurity poten- 
tial / (onsite-level difference between impurity and Fe 
atoms) and occupancy at the impurity site nimp for the 
three impurity concentrations. We found that the / val- 
ues have weak concentration dependence and were es- 
timated to be —0.4 eV for Co and —0.9 eV for Ni. 
The occupancy nimp is nearly 7.0 for Co, and 7.7 for Ni, 
apparently suggesting that extra d electrons are trapped 
at the impurity site. This is consistent with the result in 
Ref..^^ However, real-space analysis is not convenient to 
see the carrier doping effect. By visualizing the evolution 
of the Fermi-surface volume, below we will estimate the 
amount of carriers doped into the system effectively. 



Table I. Estimated onsite impurity potential I (eV) and occu- 
pancy of the impurity site Wimp for LaFei—jjMjjAsO. 



the cell is 18, so that we have 90 d bands. The panels 
(b) and (d) show spectral functions A{'k,u!) [Eq. (1)] for 
the normal cell obtained after the BZ unfolding. The ef- 
fect of impurity doping appears as small gaps or shadow 
bands in the band dispersion. As we will see later, the 
position of Ep also shifts. For the same doping concen- 
tration (5.6%), the effect of Ni-doping is more significant 
than that of Co-doping. However, for both the systems, 
the band structure does not change drastically. This is 
because the size of the impurity potential (< 1 eV) is 
still small compared to the bandwidth of the d bands 
4.5 eV). 

We show in Fig. 2 electronic density of states p{uj) — 
J2kj ^{^ — Ejk) for LaFei_2:C0:rAs0 (upper panel) and 
LaFei-ajNij^AsO (lower panel) with various impurity con- 
centrations. From the figure, we see that p{uj) around 
Ep (—0.5 eV to 0.5 eV) shows a monotonic downward 
shift with increasing impurity concentration. This behav- 
ior can be understood in terms of the center of mass of 
p{uj); since the onsite potential of Co and Ni are deeper 
than that of Fe, the center of p{uj) for the impurity sys- 
tem shifts downward compared to the undoped system. 
As the impurity concentration increases, this shift be- 
comes larger. We note that the concentration dependence 
of p{bj) is more significant for Ni than Co, due to the fact 
that the former onsite potential is deeper than the latter 
one (see Table I). Indeed, p{ijj) of the Ni-doped system 
is strongly distorted from that of the pure system. 
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Recently, Ku et al. presented a scheme to unfold bands 
in the first BZ for the supercell by utilizing the Wannier 
functions.^' The spectral function is calculated as 

A(k,c^) - -i VlmG„„(k,L.), (1) 

TT — ' 

n 

where k is wave vector in the unfolded BZ, uj is fre- 
quency, and G'„„(k, uj) is the Green function in the Wan- 
nier basis n in the primitive cell. Unless the impurity 
changes the electronic structure drastically, the quantity 
ImG„„(k, uj) can be obtained from the supercell calcula- 
tion including the impurity as 

ImG„„(k,w) = - Ejk), (2) 

KJ 

where | JK) is the J-th Bloch state with wavevector K 
in the supercell BZ and -Ejk is the corresponding eigen- 
value. The expression for (rtk| JK) is found in Rcf..^^ 

Now, let us look into the results for impurity doping. 
Figures 1 (a) and (b) [(c) and (d)] describe the Co (Ni) 
case. In the panels (a) and (c), the band dispersions for 
the 3x3x1 supercell (corresponding to 5.6 % doping) are 
plotted. The total number of transition-metal atoms in 




a)(eV) 

Fig. 2. (Color online) x dependence of the DOS p{(^) in units 
of (States/eV): (Upper panel) LaFei_j;Coa;AsO; (lower panel) 
LaFei_2;Ni2,AsO. The base lines for x = 0.056 and 0.125 are 
shifted by 40 and 20, respectively. 



In Fig. 3, we show how the Fermi surface changes by 
Co doping (left two panels) and Ni one (right two pan- 
els). For 5.6, 12.5 and 25 % concentrations, we compare 
A(k, bj ~ Ep) obtained by the BZ-unfolding method (left 
side in the two panels) with the Fermi surface obtained 
by the shift of the Fermi level of the pure system (right 
side). We see that the size and shape of Fermi surfaces 
are quite similar to each other. This clearly indicates 
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LaFeo.944Coo.o56AsO 
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(b) unfolded 
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Fig. 1. (Color online) (a) Band structure of the effective model for the d bands of LaFeo.944Coo.056 AsO in a 3x3x1 supercell. The band 
dispersion reproduces that of GGA. Due to band-folding, there are 90 bands in the first BZ. (b) Unfolded band structure. Due to the 
impurity effect, there are many gaps and several shadow bands in the dispersion, (c) A plot similar to panel (a) for LaFeo.944Nio.o56AsO. 
(d) Unfolded band structure. The size of gaps and the number of shadow bands are larger than those in the Co case. 
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Fig. 3. (Color online) Fermi surface evolution by impurity-doping obtained by the BZ unfolding method (Co: loft two panels; Ni: right 
two panels). The figures compare A(k, i?i?) in Eq. (1) with the Fermi-surface obtained by the rigid-band shift approximation. Doping 
concentration is 5.6, 12.5 and 25 % from top to bottom. 
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that substitution of Co/Ni impurities dope carriers into 
the system effectively. We can expect that such Fermi- 
surface evolution exists not only in Lallll system, but 
also in other systems such as Bal22. 

Finally, let us discuss the impurity effect on supercon- 
ductivity. Among various pairing symmetries, the pos- 
sibility of pairing with the aig symmetry has been ex- 
tensively studied. When the onsite pairing mediated by 
orbital/charge fluctuations is dominant, the gap function 
does not have sign changes in momentum space (the so- 
called s++-pairing).^^' On the other hand, the gap func- 
tion has sign changes (the so-called s±-wave) when the 
off-site pairing mediated by spin fluctuations is domi- 
nant.16) 

Experimentally,^' as mentioned above, Tc seems to 
depend only on the number of extra electrons of the tran- 
sition metal impurities. In the s++-pairing scenario, due 
to the Anderson theorem, the pair breaking effect by non- 
magnetic impurities is negligibly small. Therefore, if 
we assume that the s++-pairing interaction is suppressed 
by the Fermi surface evolution, we can naturally under- 
stand why the reduction of Tc depends only on the num- 
ber of extra electrons. 

On the other hand, the Fermi-surface evolution can be 
fatal also for the s± pairing, especially when the Fermi 
surface deforms too much and spin fluctuations are sup- 
pressed. Therefore, to explain the experiments, the pair 
breaking effect for the s-t-pairing should be also weak. 
Regarding this problem, Ref.^''^ has shown that the pair 
breaking effect is appreciable for the s± pairing. How- 
ever, the following points should be noted: First, the pair- 
breaking effect depends sensitively not only on the depth 
of the scattering potential, but also on its sign.^''^ The 
superconductivity is very fragile against impurity dop- 
ing for positive scattering potential, but rather robust 
for negative potential. As seen in Table I, Co- and Ni- 
impurity potential are negative. We should also note that 
the competition between the on-site and off-site pairing 
depends sensitively on parameters/terms in the effective 
model and how wc solve that model. In the random-phase 
approximation (RPA) for the Hubbard model, the con- 
tribution of the on-site component is estimated to be 
small. "'^^^ However, with treating the electron correlations 
beyond the RPA level or introducing electron-phonon 
couplings in the model, ^'^^ we may have situations where 
the on-site contribution to the s±-wave pairing becomes 
larger. It is an interesting future problem to examine 
these points quantitatively. 

In summary, we have studied how Co- and Ni- 
substitution change the electronic band structure of the 
iron-based superconductor, LaFeAsO. We found that 
mobile carriers are indeed doped into the system, while 
the charge density distribution seemingly suggests that 
extra electrons are trapped around the impurity atoms. 



We demonstrated unambiguously such effective carrier 
doping by visualizing the Fermi-surface evolution. The 
doping effect is well described by a rigid-band shift of 
the Fermi level. Based on these results we discussed how 
the s-wave (s±-wave) scenario can be compatible with 
experiments. 
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